Correlated electron transport through parallel double-quantum-dot 



IT) 

o 
o 

(N 
5h . 

< 
(N 



13 

o 
o 



> 
oo 
oo 
(N 

o 
in 
o 

-I— > 



o 
o 



X 



Rong Lii 1 , Zhi-Rong Liu 2 , and Guang-Ming Zhang 2,1 

1 Center for Advanced Study, Tsinghua University, Beijing 100084, China; 
2 Department of Physics, Tsinghua University, Beijing 100084, China 
(Dated: February 2, 2008) 

We investigate the spectral and transport properties of parallel double-quantum-dot (DQD) sys- 
tem with interdot tunneling coupling in both the equilibrium and nonequilibrium cases. The special 
geometry of DQD system is considered, in which each dot is connected to two leads by the tunneling 
barriers. With the help of Keldysh nonequilibrium Green function technique and the equation-of- 
motion approach, the spectral function and the conductance spectra of DQD system are calculated 
in two cases with and without the intradot Coulomb interaction, respectively. The exact calculation 
is performed in the absence of intradot Coulomb interaction. For the case with intradot Coulomb 
interaction, the Hartree-Fock approximation is applied to truncate the equation of motion for the 
high-order Green functions at high temperatures. The phenomenon of correlated electron transport 
is clearly shown in the linear conductance of each dot in the presence of interdot tunneling when 
setting one dot level and tuning another. The interplay between the intradot Coulomb interaction 
and the interdot tunneling coupling is displayed. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73.21. La 



I. INTRODUCTION 

Recently, due to the rapid progress in nanotech- 
nologies, quantum transport through DQD system has 
been the subject of active theoretical and experimental 
research, 1 Compared to the single quantum dot, the in- 
terdot coupling and the intradot on-site Coulomb inter- 
action in DQD system could generate novel many-body 
states, e. g. the molecular Kondo stated Apart from 
its importance in understanding some basic problems in 
condensed matter physics, the DQD devices are also cru- 
cial ingredients in the emerging field of spintronics and 
several quantum computation schemes designed with the 
electron spins or with the coherent mode in an artificial 
moleculei^ 

Transport through DQD system in the Coulomb 
blockade and the Kondo regimes has already received 
some theoretical attention, in which two dots are ar- 
ranged in series, 5 ! 6 ! 7 ! 8 ! 9 ! 10 ' 11 ! 12 ' 13 ! 14 ! 15 ! 16 paralleli&iii 8 . or 
T-shapeii2i22i2i In most of these geometries, the conduc- 
tance spectra have been studied for electron transport 
through the whole DQD system, while the electron trans- 
port through each individual dot could not be detected. 
In this paper we consider a new configuration of DQD 
system (see Fig. 0, in which each quantum dot is con- 
nected to two leads by the tunneling barriers. In this 
geometry, the collected electron transport phenomenon 
can be investigated by measuring the transport proper- 
ties through one individual dot when tuning the level of 
another dot and its interdot coupling. The similar con- 
figuration of DQD system has been proposed in Ref. 12^ . 
however, the detailed investigation on the spectral and 
conductance spectra, and the interplay between interdot 
tunneling coupling and intradot on-site Coulomb inter- 
action has been not discussed yet. 

In this paper, by using the Keldysh nonequilib- 
rium Green function technique and equation-of-motion 



method, we calculate the spectral and the conductance 
spectra of tunneling-coupled DQD with and without in- 
tradot Coulomb interaction, respectively, in both the 
equilibrium limit and in the nonequilibrium case. In Sec. 
II, we formulate the model Hamiltonian of the special ge- 
ometry of DQD system, and deduce the current formula 
through each quantum dot. 

In Sec. Ill, for DQD without the intradot Coulomb in- 
teraction, we exactly calculate the Green functions, from 
which both the spectral function and the conductance 
are obtained, showing the resonant peaks corresponding 
to the molecular resonant-tunneling bonding and anti- 
bonding states. The image of total differential conduc- 
tance shows " anticrossing" phenomenon when the lev- 
els of two dots match. It is interesting to find that in 
the presence of interdot tunneling coupling, the electron 
transport through one dot is strongly influenced by the 
electron transport through another dot, namely, when 
setting dot 1 level (ei) and tuning dot 2 level (£2) the 
linear conductance of dot 1 (Gdoti) shows a dip at the 
same position of £2, where the linear conductance of dot 
2 (Gdot2) displays a resonant peak. Increasing the inter- 
dot tunneling coupling gives rise to a more complicate 
structure of linear conductances through each dot. 

In Sec. IV, we apply the Hartree-Fock approximation 
to DQD with the intradot Coulomb interaction to trun- 
cate the equation of motion for the high-order Green 
functions, which is known to capture the correct qual- 
itative feature of physics of single quantum dot in the 
Coulomb blockade regime. Under the same approxima- 
tion, we also obtain the Keldysh less Green function, 
which is needed in the self-consistent evaluation of the 
dot occupation numbers when a finite bias voltage is ap- 
plied across the system. As the intradot Coulomb in- 
teraction becomes the largest energy scale in the single- 
electron tunneling regime of quantum dot, the intradot 
Coulomb interaction leads to two groups of peaks in 
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both the spectral and conductance spectra, and the in- 
tcrdot tunneling coupling causes the splitting of peaks in 
each group. Compared to the case without the intradot 
Coulomb interaction, when setting e\ and tuning £2, one 
more dip in Gdoti and resonant peak in Gdoti will ap- 
pear corresponding to the Coulomb interaction energy, 
which also provides the clear evidence of the correlated 
electron transport through DQD system. The paper is 
closed with a brief summary in Sec. V. 



II. PHYSICAL MODEL AND CURRENT 
FORMULA 



dot 1 




dot 2 



Quantum dots behave as artificial atoms and single- 
electron transistors in their charge and energy quantiza- 
tions due to the small dimensions (compared to the Fermi 
wavelength). They are often described by the single- 
impurity Anderson model, in which there is a Coulomb 
repulsion between electrons in the dot. The parallel DQD 
system could be modelled by using a two-impurity Ander- 
son model with an extra interdot tunneling corresponding 
to the electrons hopping between two dots. Each QD is 
attached to two electronic leads with different chemical 
potentials by quantum-tunneling barriers. Geometry of 
the system is shown schematically in Fig. ^ The model 
Hamiltonian is then given by 



H = H, 
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DD 
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(1) 



The first term Hi eac i s describes the noninteracting elec- 
trons in the 77-th lead with the electron creation and 
annihilation c^, ia operators, 



k,?7,o" 



(2) 



where 77 = 1, 2, 3, 4 corresponding to the four leads shown 
in Fig. n k denotes the wave vector, and a is the spin 
index. The second term Hdd in Eq. (JIJ describes the 
electrons in the dots 1 and 2 with the interdot tunneling 
coupling t and the intradot on-site Coulomb interaction 
U, 



Hdd — 



U ■ 



(■2d 2a d 2 <, + — > , n 2(7 n-m 
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(3) 



The last term in Eq. Q describes the tunneling coupling 
between dots and the external leads with the tunneling 
matrix elements V, 



k , <j 

+ E f^cL^dto + V A cl 4a d 2a + h.c.) . (4) 



FIG. 1: Schematic diagram of the parallel DQD system stud- 
ied in this paper. 



When we use the Keldysh nonequilibrium Green func- 
tion technique to express the current through dots 1 and 
2, the respective conductance can be written in terms of 
the distribution functions of leads and the local proper- 
ties of DQD system. To this end we write the current 
from the lead 1 to dot 1 as 



Ji 



-e (dNx/dt) = - (ie/h) ([H, Nx}) 



(5) 



where N% = J^w a c ki CT c ki<r- With the help of Eq. QJ, 
one can easily obtain the current from the lead 1 to dot 
1 as 



where the Keldysh less Green functions are defined as 
Gf aMa (t,t') ee i(ct la (t')d la (t)), 



(7) 



By using the Keldysh Green function formalism, one 
can obtain the Dyson equations as 



~9kla M G \*,l<r ( W )] , 
-glla M G< ( W ) 



(8) 



where the Keldysh less Green function for electrons in 
dot 1 is defined as 



k.cr 



G< M (t,t') = i(dl {t')d la (t)), 



(9) 
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and (?kio- is the Green functions for noninteracting elec- 
trons of the lead 1. The time-ordered and anti-time- 
ordered Green functions are denoted by the superscripts 
t and t, respectively. Therefore, the current from the lead 
1 (or 2) to the dot 1 can be expressed as 

+/i(2) (e) [G[ aM (e) - G\ aAa (e)] } , (10) 

where the line-width I\(2) — ^PfVim) with pp the den- 
sity of states of leads, and fi(2) ( e ) is the Fermi-Dirac 
distribution functions of the lead 1 (or 2) with the chem- 
ical potential p\{2)- Here we have assumed that the leads 
give rise to a fiat, energy independent, density of states 
(i. e., the wide-band limit). By applying the similar pro- 
cedure, the currents from the leads 3 and 4 to the dot 2 
are found similarly, 

J 3( 4 ) = f E/ ^ r 3(4) (e) {G^ 2(T (e) 

(J 

+/ 3 (4) (e) [G r 2aM (e) - G a 2a , 2a (e)] } , (11) 



with r 



3(4) 



for the lead 3 (or 4). 



Under the symmetric condition: \i\ + /i 2 = 1^3 + Hi, 
there is no current between dot 1 and dot 2. In the 
steady state, the current will be uniform, so that the 
current through dot 1 satisfies Jdoti = Ji = —J2, and 
then one can symmetrize the current through the dot 1 
as 



Jdo 



J2) 



= ^E/ ^[(f^~m(G r laM (e) 

(7 

-Gl aila {e))+^i-Ti)Gf aM {e)], (12) 
and similarly the current through the dot 2 is 

J dot2 = ^E/ ^[(/3r3-.Ar 4 )(G^ 2CT ( e ) 

£7 

-GL,2 ff (<0) + (T 3 - T 4 ) G 2 < >2CT (e)] . (13) 

It is noted that Eqs. (|12f> and l|13|) give rise to the cur- 
rent through each dot of the DQD system in terms of the 
distribution functions of leads and the local properties of 
dots. In order to obtain the current, one has to compute 
the retarded and the Keldysh less Green functions of the 
DQD system in the presence of both the interdot tun- 
neling coupling t and intradot Coulomb interaction U as 
well as the tunneling coupling of DQD system into the 
leads. 

Without the loss of generality, we assume that the dot- 
lead tunneling coupling are symmetric: ri = T 2 = T3 = 
T4 = r, and the symmetric configuration of the chemical 
potentials: (Xi = /Z3 = fi+eV/2, and p 2 = Hi = fi—eV/2. 



Then the current through the dot 1 (or 2) becomes 
Jaotm = § E J de [h (e) - h (e)] A 1[2)a (e) , (14) 

IT 

where the dot spectral functions are defined as 



A 1{2)a = (-l/ir)Im G r 1(2)tjM2)a 



(15) 



The associated differential conductance at /1 = is found 
to be 

dJ do ti(2) e 2 (3T ^ f 

CJ 

(16) 



(A<-*) + i) (e^ + *) + l) 

and the linear conductance at zero bias is thus given by 

r e2T sr I a ( dfFD ^ 

G dotl{2) = —}^ J de \- 



ck 



A m *{e), (17) 



where fpu (e) is the Fermi-Dirac distribution function. 



III. THE NONINTERACTING DOTS 

In this section we study the case of noninteracting dots, 
i. e., U = 0. In this case, we can exactly derive the 
Green functions, the associated spectral, and transport 
properties of the DQD system, which shows the clear 
evidence of correlated electron transport through DQD. 

By applying the equation-of-motion approach, we have 



(iu n - e x ) {(d 1(T \d\ c 
l-t((d 2a \dl))+V 1 *J2(( c ^\ d i 



(18) 



- £ kl) ((c k l ff |4a)) = Vl ((dl ff |4a)) . ( 19 ) 



(icj n - e k2 ) ((ck2a\d\ a )) = V 2 {{d la \d\ 



(20) 



Here we assume that the conduction leads have a flat and 
energy independent density of states (i. e., the wide-band 
limit), leading to 
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FIG. 2: The images of total linear conductance Gdoti + Gdot2 
as a function of the dot levels ei and ti for noninteracting 
(U = 0) double dots with V = 0.01, k B T = 0.01, t = in (a), 
and t = 0.1 in (b). 



and then 



iu>„ ~e 1 + -(T 1 + T 2 ) 



dia\d[ c 



t((d 2 a\d\ c 



1. (21) 



By applying the similar method, we obtain the equation 
of motion for the Green function ( (e^o-Mio-// as 



e 2 + 2 ( r 3 + r 4 ) 



From the above two equations and performing the an- 
alytic continuations, we can deduce the retarded Green 
functions in the symmetric dot-lead tunneling coupling 

case, 



d2*\d\ a/ 
t((dia\dl)) = 0. (22) 
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FIG. 3: The spectral function as a function of frequency in 
(a), and the differential conductance as a function of bias 
voltage in (b) with the dot-level: ei = £2 = for different 
interdot tunneling couplings. 
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Gr 
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where 



£±= 2 



(ei + e 2 ) ± V (ei - e 2 



4i 2 



(23) 



(24) 



denote the energies of antibonding and bonding molec- 
ular resonant-tunneling states, respectively. Note that 
the bonding state moves down in energy with the inter- 
dot tunneling coupling and the antibonding state moves 
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up. The energy difference between the antibonding and 

bonding states is e + — e_ = J (ei — £2) 2 + 4i 2 . When 
the levels of two dots cross, i. e. t\ = £2, one has an an- 
ticrossing of e+ and e_ with the minimum antibonding- 
bonding energy difference 2t. For large dot-level energy 
difference, the eigenenergies of tunneling-couplcd double 
dots approach to the energy levels of uncoupled dots, £1 
and 62. 

From the retarded Green functions, the spectral func- 
tion of dot 1 is found to be 



Ala H 



(2 



r 2 



(-2 



J (we-)' 



r 2 



(25) 



and changing the index 1 <-> 2 gives the spectral function 
of dot 2. The spectral functions of dot 1 and 2 display two 
Lorentzian resonances at the same positions: e±, while 
the ratio of their heights is (e+ — £2) / (ea — £-) for the 
dot 1, and (e + — ei) / (ei — e_) for the dot 2. 

In Fig. |2| we show the images of the total linear conduc- 
tance Gdoti+Gdot2, obtained by Eq. JHJl, versus changes 
of levels ei and £2 of two dots. The bright regions corre- 
spond to high conductance and the dark regions to low 
conductance. For completely decoupled dots [t = 0) as in 
Fig. |2[a) , tuning the level of one dot changes the charge 
on this dot without affecting the charge on the other dot. 
For nonzero interdot tunneling coupling, a similar "anti- 
crossing" occurs at t\ = £2 in the diagram of conductance 
[see Fig. Hb))]. 

In Fig. |21 the spectral function and the differential 
conductance of dots are displayed for e% = £2 = and 
different interdot tunneling couplings. In the absence of 
interdot tunneling, the electron transports through each 
dot independently, which exhibits a resonant peak when 
the dot levels match the chemical potential of the leads. 
Increasing the interdot tunneling leads to the splitting of 
resonant peak at the positions e±, corresponding to the 
bonding and antibonding molecular states of the DQD 
system. Tuning the bias voltage gives rise to the reso- 
nances in the differential conductance when the resonant 
peaks in the spectral function enter the region between 
the chemical potentials of leads. 

We also show the spectral function and the differential 
conductance of dots along the direction t\ = £2 for differ- 
ent interdot tunneling couplings in Fig. 0] Compared to 
Fig. |3 the spectral function is shifted to ei — £2 =0.1. 
The asymmetric spectral function along to leads to the 
four peaks structure in the differential conductance as a 
function of the bias voltage for the general configuration 
of dot level and the finite interdot tunneling coupling [see 
the green line of Fig. E£b)]. One interesting observation 
is that under the condition t — t\ — £2, this general four 
peaks structure reduces to the three peaks structure [see 
the red line of Fig. Efb)] > while in this case the antibond- 
ing energy shifts to zero [see the red line of correspond- 
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FIG. 4: The spectral function (a) and the differential con- 
ductance (b) with the dot-level: ei = £2 =0.1 for different 
interdot tunneling coupling. 



ing spectral function in Fig. Ufa)]. This feature offers a 
possible method to measure the tunneling coupling be- 
tween dots in experiments. Fig. [5] displays the spectral 
function and the differential conductance along the di- 
rection t\ = — £ 2 = — e c i for different interdot tunneling 
couplings. In this case, from Eq. I|24|) one can easily 
see that the bonding and antibonding energies reduce to 
e± = ± (td + 1), which yields two resonant peaks in the 
spectral function for finite interdot tunneling coupling, 
while the heights of these two peaks are different for the 
ratio of heights is t/ (2ed + t) for dot 1 and (2ed + t) /t 
for dot 2. As these two resonant peaks sit at the posi- 
tions with the same absolute value, there are only two 
resonances in the differential conductance compared to 
the four resonances in the case £1 = £2 (see Fig. |SJ). 

In Fig. El the linear conductance of each dot as a 
function of one dot's level is displayed by setting an- 
other dot's level. In the absence of the interdot tunneling 
coupling, electron transports independently through each 
dot, which shows one constant linear conductance of dot 
1 and one resonant peak in the linear conductance of dot 
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FIG. 5: The spectral function of the dot 1 in (a) and the 
dot 2 in (b) and the differential conductance in (c) with the 
dot-level: ei = —0.1,62 = 0.1 for different interdot tunneling 
coupling. 



FIG. 7: The linear conductance of dot 1 (a) and 2 (b) as a 
function of dot level of one dot 62 by setting another dot's level 
as ei = when increasing the interdot tunneling coupling, and 
the associated spectral functions for different 62 [(c) and (d)]. 
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FIG. 6: The linear conductance of dot 1 and 2 as a function of 
dot level of one dot 62 by setting another dot's level as ei = 
without (a) and with (b) the interdot tunneling coupling. 



2 when £2 matches the chemical potential of leads. In the 
presence of interdot tunneling coupling, the linear con- 
ductance of dot 1 shows a dip at the same position of £2, 
where the linear conductance of dot 2 exhibits a resonant 
peak. 

An interesting feature is shown in Fig. [7|when increas- 
ing the interdot tunneling coupling: the resonant peak in 
the linear conductance of dot 2 is suppressed into two 
splitting peaks. Both the dip in Gdoti and the split- 
ting peaks in Gdot2 could be understood with the help 
of the associated spectral function of each dot for differ- 
ent dot level €2- The blue line in Au2)a- °f Fig- s h° ws 
—dfF d /de in the formula of linear conductance Eq. I|17(l , 
which sets the window of temperature. It is noted that 
the spectral functions of both dots will include two peaks, 
corresponding to the bonding and antibonding states at 
finite interdot tunneling coupling, while the heights of 
these two peaks are different. For the spectral function 
of dot 1 {Ai a ), when increasing the dot level €2, the peak 
corresponding to bonding state e_, which is smaller than 
the minimum of ei and 62 , moves to the window of tem- 
perature with the increasing of its height. This will en- 
hance the linear conductance of dot 1 as increasing £2 
from 62 = 0, and the linear conductance saturates to a 
constant value as the bonding-state peak completely en- 
ter into the window of temperature [see Fig. EJ a )]- For 
the spectral function of dot 2 {A<2 a ), as increasing the dot 
level 62, the bonding-state peak moves into the window 
of temperature, but the height of this peak decreases. 
This leads to a maximum in linear conductance of dot 
2 [see Fig. Efb)]. While this maximum in the linear 
conductance of dot 2 could not be observed in the small 
interdot coupling case (as shown in Fig. [|J| or in the 
low temperature case, because the bonding-state peak is 
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too close to lj = or the window of temperature is too 
narrow. 



IV. THE INTERACTING DOTS 

In this section we study the transport phenomena 
through tunncling-coupled DQD system with intradot 
Coulomb interaction with the help of equation-of-motion 
approach. This method consists of differentiating the 
Green function with respect to time, thereby generating 
higher-order Green functions in the presence of intradot 
on-site Coulomb interaction which eventually have to be 



closed. Here we apply the Hartree-Fock approximation 
to truncate the higher-order Green functions, which in- 
cludes the contributions of higher order electron tunnel- 
ing processes. It is known that the Hartree-Fock approxi- 
mation captures the correct qualitative feature of physics 
of single quantum dot in the Coulomb blockade regime, 
valid in the higher temperature regime. 

By applying the equation-of-motion approach and the 
Hartree-Fock approximation to truncate the higher-order 
Green functions (for detailed calculation see the Ap- 
pendix A), we obtain the Green functions for dot elec- 
trons as 



fi [i 



iui n 



£2 



(r 3 + r 4 



lUJ n - 61 



lUJ n — 62 + t 2 J 



t 2 F 1 F 2 



(26) 



d 2a \d\ a 



F 2 [iuj n -e 1 + F 1 i(T 1 +T 2 ) 



, jp i(ri+r 2 ) 



iuj n — e 2 + t 2 ^ 



t 2 F x F 2 



(27) 



and 



d 2a \d\ a )) = Ud la \d\„ 



-tF\F 2 



, jp »(ri+r 2 ) 
tuj n -ei + Fx^-^ — '- 



lLU n — 6 2 + t 2 2 



(28) 



where 

Fi(uj n ,ei,{nm)) = 1 + - =, 

lUJ n - 61- U 

F 2 (1^,62,(71^)) = 1+ . UM - ■ (29) 

ZW„ —6 2 — U 

Performing the analytic continuations, we can deduce 
the retarded Green functions for dot electrons in the 
tunncling-coupled DQD system with intradot on-site 
Coulomb interaction: G\ a j a , where i,j — 1,2 corre- 
sponding dot 1 and 2 respectively. The occupation num- 
bers are subjected to the self-consistency condition: 

{***} = -i J f^G&^M, (30) 

where i = 1 and 2. 

Next we calculate the Keldysh less Green functions 
Gf aia (cl>), which is needed in the self-consistent equa- 
tion Eq. (|3UII for the dot occupation numbers. Under 
the same Hartree-Fock approximation as in the retarded 
Green functions, the equation of motion approach gives 
that the Keldysh less Green functions satisfy (for detailed 
calculation see the Appendix B) 

Gf aM = i[r 1 / 1 ( w ) + r 2 / 2 ( w )]GL )lff G? CTilCT 

+l [r 3 / 3 M + r 4 / 4 H] G r 2aM G a 2aM m 



for dot 1, and 

G 2cr , 2a = i[Tifi(Lj)+T 2 f 2 (c)]Gl a<2a G a i aj2a 

+l [r 3 / 3 M + r 4 / 4 (w)] GL, 2CT GL, 2ff ,(32) 

for dot 2. 

Fig. IHJshows the images of the total linear conductance 
Gdoti + Gdot2 of DQD with the intradot on-site Coulomb 
interaction versus the changes of the dot-level 6\ and e 2 . 
Without interdot tunneling coupling, the images of con- 
ductance shows a square lattice structure [see Fig. ISta)]. 
In the presence of both the interdot tunneling coupling 
and the intradot Coulomb interaction, the " anticrossing" 
occurs at the positions of ei, 6 2 , a + U, and 6 2 + U [see 
Fig. Eft)]. 

The effects of the intradot Coulomb interaction and 
the interdot tunneling coupling can be clearly shown in 
the spectral functions of quantum dots. Fig. |5| is the 
spectral function and the differential conductance of the 
dot 1. In the absence of the interdot tunneling coupling, 
the intradot Coulomb interaction leads to two peaks at 
£i(2) and ei(2) + U in the spectral function, in agree- 
ment with the well-known result of Coulomb blockade at 
higher temperatures. These two peaks split in the pres- 
ence of the interdot tunneling coupling, and the spectral 
function displays a four peak structure. The differential 
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FIG. 8: The images of total linear conductance Gdoti + Gdoti 
as a function of the dot levels ei and £2 for interacting (U = 1) 
double dots with T = 0.01, k B T = 0.01 t = in (a), and 
t = 0.1 in (b). 



conductance exhibits the resonance peaks when the spec- 
tral function peaks enter the window of bias voltage. Due 
to the small size of QD, the Coulomb charging energy U 
becomes the largest energy scale. In the case U > t, the 
peaks in both the spectral function and the linear con- 
ductance form two distinct groups separated by roughly 
the intradot Coulomb interaction U, while in each group 
there are two peaks separated by roughly the interdot 
tunneling coupling t. 

In Fig. we show the linear conductance of each 
dot with the intradot Coulomb interaction as a function 
of one dot-level by setting another dot-level. In the ab- 
sence of interdot tunneling coupling, by setting ej and 
tuning £2, Gdoti remains constant and Gdoti shows two 
resonant peaks when the levels £2 and £2 + U match the 
chemical potential of leads. Due to the interdot tunnel- 
ing coupling, the electron transport through one dot is 
strongly influenced by the electron transport through an- 
other dot: Gdoti displays two dips at the corresponding 
resonant positions of Gdoti- This feature provides a pos- 
sible experimental method to observe the correlated elec- 
tron transport phenomena in the real interacting DQD 
system with strong intradot tunneling coupling at high 
temperatures. 



> 

T3 




FIG. 9: The spectral function (a) and the differential con- 
ductance (b) of the dot 1 in the DQD system with the in- 
tradot Coulomb interaction in the general case of the dot- 
level: £1 = —0.1, £2=0 for different interdot tunneling cou- 
pling. 



V. SUMMARY 

In this work we have investigated theoretically the elec- 
tron transport through parallel tunneling-couplcd DQD 
system with the special configuration that two leads are 
connected to each dot. By applying the Keldysh nonequi- 
librium Green function technique and the equation-of- 
motion approach, we have obtained the spectral and con- 
ductance spectra of DQD with and without the intradot 
on-site Coulomb interaction, respectively. 

The main results are summarized as following: 

1) In the absence of intradot Coulomb interaction, the 
exact results of the spectral function and the conductance 
and obtained, showing resonant peaks at the positions 
of antibonding and bonding molecular states for finite 
interdot tunneling couplings, while the heights of these 
resonant peaks depend on the levels of two dots. 

2) Tuning the dot levels results in different structure 
of the differential conductance, e. g., under the condition 
£1 = £2 = t, the general four-peak structure in the dif- 
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FIG. 10: The linear conductance ol dot 1 and 2 with the 
intradot-dot Coulomb interaction as a function of dot level of 
one dot £2 by setting another dot's level as ei = without [in 
(a)] and with [in (b)] the interdot tunneling coupling. 



ferential conductance can reduce to the three-peak struc- 
ture. 

3) For DQD with the intradot Coulomb interaction, 
the Hartree-Fock approximation has been applied to cal- 
culate both the retarded and the Keldysh less Green func- 
tions, which gives the transport properties at high tem- 
peratures. 

4) As the interdot Coulomb interaction U becomes the 
largest energy scale of quantum dot in the single-electron 
tunneling regime, the peaks in the spectral function and 
the conductance form two groups separated by roughly 
{/, and the interdot tunneling coupling t induces the split- 
ting of peaks in each group. 

5) The linear conductance of each dot shows the 
phenomena of correlated electron transport through 
tunncling-coupled DQD system, which could be observed 
with the present experimental techniques. 
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APPENDIX A: THE EVALUATION OF 
RETARDED GREEN FUNCTIONS 

In this Appendix, we give the detailed evaluation of 
the retarded Green functions for dot electrons in the 
tunncling-coupled DQD with intradot on-site Coulomb 
interaction based on the equation-of-motion approach. 
Its equation of motion is yielded as 



(icu„ - ei) Ud 1(7 \dl c 
= 1 - t ((d 2ry \d\ a )) + U ((n lw d 1(7 \dl 



(iuj n - e k i) Uckiv\d\ a )) = Vl \(d llT \d\ c 
(iuj n - e k 2) ((ck2ffK<r// = V 2 {{d la \d\ c 



(A2) 



(A3) 



where the higher-order Green function (J^ni-^di^d^ 

is generated. The equation of motion of (^nivdi,j \d\ a 
can be further deduced to 

{iu n -ei-U) ^n w d la \d\ a )) 
= {nia} - t r/ZniF^o-lrfitr 
- ((dU lw d lr7 \dl))) 
"E( y i (( c lmdmd la \d\ 

k 

-V* ((d\ w c km d la \d\ 



d\ w d 2 adia\d\ a 

V* ((niaCkXa\d\ 



(A4) 



Making the Hartree-Fock approximation to decouple the 
higher-order Green functions in the above equation: 



la), 


> - (4^w)< 


X^Miff); 


>, 
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> - (4wdm)( 


;(^i4ff); 
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L), 




((rfiffML) 


>• 


iff). 


) ^ (d[ w c km ) 


((dlMa) 


)(A5) 



10 



and considering the facts: Vi = V{ , (d\—d2w 
d^dm), and (c^dw) = /d\ w c kl a), the higher-order 



Green function is thus derived as 

[ni w dia\d\ t 
(niw) 



ei - u y 



l-t((d 2 a\d\ c 



E^^ + E 

Introducing 

Fi (iw„,ei, (nia)) = 1 + 

F2 («W n , 62, (n2a)) = 1 + 



4W„ - £kl ^ l^n - ek2 



IV2 



(A6) 



^ (nig) 
iw„ — ei — U' 

U (n 2 g) 
iuj n - e 2 - U ' 



(A7) 



the dot Green function ((di a \d\ a JJ is found to satisfy 
the following equation: 



IV2 



IW„ - €kl l^n - £k2 



+tF 1 ((d 2a \dl)) = Fx. (A8) 



By applying the similar method, the equation of motion 
of the Green function \(d2a\d\ (T )) satisfies: 



tF 2 ((d la \dl 



Jn - e 2 - F 2 E 



E 



4W„ — £k3 ^ iW n — ek4 



0. 



(A9) 



From the above two equations, we obtain the retarded 
Green functions as shown in Eqs. H2tjf> an d (EH- 
Similarly, we can obtain the retarded Green functions 

(^c^o-ldL)) and ^^laML)) as shown in Eqs. 

and (gill- 



APPENDIX B: THE EVALUATION OF THE 
KELDYSH LESS GREEN FUNCTIONS 



needed in the self-consistent evaluation of the dot oc- 
cupation numbers when a finite bias voltage is applied 
across the system. Its equation of motion is yielded as 



= 61* (4„ (0 Ala (t)) ~ ti (d\ a (f) d 2a (t)) 

+Ui (d[ a (<') [n m d la ] (t)) + ^1* E 1 (0 c ki CT (*) 

k 

(Bl) 



and then we perform the Fourier transform 



vGfv ,1a M 



- ^G< aM ( W ) - tG&. lff (a;) + U 2 * £ G< ^ (0 

k 

-c/ (UMdi^Y M + E G ki.,i^ M ( B2 ) 



Under the same Hartree-Fock approximation as in the 
retarded Green functions, the equation of motion for the 

higher-order Keldysh Green function ((nicdia^ia 
can also be deduced to 



(w - e x - U) ((n^d la \d\ a )Y H 
= -t <n 1? ) G< aM M + V* <n w ) E G ki CT ,i. H 

k 

+^2*ME G k2M CT M- (B3) 
k 

By using the exact Dyson equation for G^^u la (to): 

G kl(2)<r,l<7 ( W ) = V l [ffkl(2)<r M H 

+.9ki (2 ).HGL,i CT H] . (B4) 



With the help of equation-of-motion approach, one can 
also obtain the Keldysh less Green functions, which is the equation of motion for the Gf a la becomes 



11 
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By applying the same method, we obtain the equation of motion for the Keldysh less Green function G^ CT lt7 as 



(B5) 



tF 2 Gf alc 



u;-e 2 ~F 2 [ \V 3 \ 2 Y,&3* 



^2<t,i<7 — -Pa 



Ga 

2fT,lc 



(B6) 



From these two equations, we can solve the Keldysh less dot 2. 
Green functions as Eq. (|31|l for dot 1 and Eq. (|32|l for 
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